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Abstract 

We show a perturbation theory of the Josephson transport through ferromagnetic insulators (FIs). Recently we have 
found that the appearance of the atomic scale 0-7T transition in such junctions based on numerical calculations. In 
order to explore the mechanism of this anomalous transition, we have analytically calculated the Josephson current 
using the tunneling Hamiltonian theory and found that the spin dependent 7r-phase shift in the FI barrier gives the 
atomic scale 0-7r transition. 
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1. Introduction 

In the usual Josephson junctions at equilibrium 
the phase difference of the superconducting order 
parameter on the two superconductor is zero. On 
the other hand, in the Josephson junctions with 
ferromagnetic- metal interlayer (S-MF-S junctions), 
the ground state may correspond to the phase dif- 
ference [1,2]. The Josephson 7r-junction formation is 
related with the damping oscillatory behavior of the 
Cooper pair wave function in a FM [3,4]. In terms 
of the Josephson relationship Ij — 1q sin 0, where 
<j> is the phase difference between the two supercon- 
ductor layers, a transition from the to n states 
implies a change in sign of Iq from positive to nega- 
tive. Experimentally the existence of the 7r-junction 
in S/FM/S systems has been confirmed by Ryan- 



zanov et al. [5] and Kontos et al. [6] Until recently, 
however, investigations on the 7r junction have been 
mainly focused on the S/FM/S systems. 

We have predicted a possibility of the 7r-junction 
formation in Josephson junctions through ferro- 
magnetic insulators (FIs) by numerically solving 
the Bogoliubov-de Genne equation [7,8] and the 
Nambu Green's function [9-12]. The formation of 
the 7T junction using such an insulating barrier is 
very promising for future qubit application [13— 
16] because of it's low decoherence nature [17,18]. 
More importantly, we have shown that the ground 
state of such junction alternates between 0- and tt- 
states when thickness of FI is increasing by a single 
atomic layer. In this paper in order to understand 
the physical mechanism of this atomic scale 0-7r 
transition, we analytically calculate the Josephson 
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current based on the tunneling Hamiltonian method 
and show that the spin-dependent 7r-phase shift of 
the tunneling matrix element of the FI layer is the 
origin of this anomalous transition. 



In this paper we focused on the high- Tc ci-wave 
junction with a FI barrier [Fig. 1(b)]. We note that 
the qualitatively same result can be obtained for the 
case of conventional s-wave junctions. 



2. Tunneling properties of a ferromagnetic 
insulator 

In this section, we briefly describe the electronic 
structure of a representative of FI materials, i.e., 
La 2 BaCu0 5 (LBCO) [19]. The typical DOS of 
LBCO for each spin direction is shown schematically 
in Fig. 1(a). The exchange splitting V cyi is estimated 
to be 0.34 eV by a first-principle band calculation 
using the spin-polarized local density approxima- 
tion [20]. Since the exchange splitting is large and 
the bands are originally half-filled, the system be- 
comes FI. Based on the band structure [Fig. 1(a)], 
the spin a(= J [, J,) dependent transmission coefficient 
T a for the FI barrier can be calculated as [7] 



Pat 
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(1) 



by using the transfer matrix method, where Pfm — 
— (+)1 and ctL F is a spin-independent complex con- 
stant, t is the hopping integral in the FI, g is the gap 
between the up- and down-spin band, and Lp is the 
layer number of the FI. We immediately find 



1) 



Lf 



(2) 



so the relative phase of T a between spin up and down 
is 7r (0) for the odd (even) number of Lp . Next sec- 
tion, we will calculate the Josephson current through 
such a FI analytically. 



(b) 
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Fig. 1. (Color online) (a) The density of states for 
each spin direction for a ferromagnetic insulator, e.g., 
LBCO, and (b) schematic picture of c-axis stack high-T c 
d-wave superconductor/LBCO/high-T c d-wave superconduc- 
tor Josephson junction. 





3. Theory 

In this section, we analytically calculate the 
Josephson current of S/FI/S junctions based on the 
tunneling Hamiltonian approach. Let us consider 
a three-dimensional S/FI/S junction as shown in 
Fig. 1(b). The Hamiltonian of the system can be 
described by 



H — Hi + H-2 + Ht + H< 



Qi 



(3) 



where H± and H2 are Hamiltonians describing the 
d-wetve superconductors: 
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where p, is the chemical potential and ip [ip ) is the 
fermion field operator and g (r — r') is the attractive 
interaction. The tunneling Hamiltonian with a spin- 
dependent tunneling matrix element t a of the FI 
barrier is given by 

H T =Y,f dr J dr ' [«<r(' , y)4w^M 

+ h.c] , (5) 

and 

(Q1-Q2) 2 



H, 



Q 



8C 



(6) 



is the charging Hamiltonian, where C is the capac- 
itance of the junction and Q\(2) is the operator for 
the charge on the superconductor 1 (2), which can 
be written as 



(T) 



By using the functional integral method [17] and 
taking into account the spin dependence of t a ex- 
plicitly, the ground partition function for the system 
can be written as follows 
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where f3 — l/fc^T, is the Grassmann field 

which corresponds to the fermionic field operator 
m Eq- (2), and the Lagrangian L is given by 

x ( r ) = EE / dr ^- ( r > r ) ( r ' r ) + 

a i=l,2 

(9) 

In order to remove the ip term in the Hamiltonian 
H(t), we will use the Hubbard-Stratonovich trans- 
formation which introduces a complex pair potential 
field 

A(r,r';r) = |A(r, r'; r)| exp [i<j> (r, r'; r)] . (10) 

The resulting action is only quadratic in the Grass- 
mann field, so that the functional integral over this 
number can readily be performed explicitly. The 
functional integral over the modulus of the pair po- 
tential field is taken by the saddle-point method. 
Then the partition function is reduced to a single 
functional integral over the phase difference <fr = 
<pi — 4>2- To second order in the tunneling matrix 
element, one finds 



= J D4>(r) 



exp 
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where the effective action is given by 
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In the calculation we have ignored the irrelevant 
quasiparticle tunneling term for simplicity. Here 



Ic 
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(13) 



is the Josephson critical current, and then the 
Josephson current is given by 



Ij(4>) = -fcsint 



(14) 



where the negative (positive) Ic corresponds to the 
7T (0) junction. The Josephson kernel j3(r) is given 
in terms of the off-diagonal Nambu Green's function 
for two superconductor Fi (i — 1, 2), i.e., 



P{t) = -t E *K fe ' fe ')*t( fe ' W ( fe ' T ) F 2 (*', ~r) 



k.k' 

The Nambu Green functions is given by 



(15) 



Fi (fc,w„) 
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where £u = h 2 k 2 /2m— is the single particle energy 
relative to the Fermi surface and fiw„ = (2n + l)-7r//3 
is the Matsubara frequency (n is an integer). In the 
case of the cuprate high-T c superconductors with 
the d x 2_ y 2 symmetry [21], the order parameter is 
given by 



Aj(fc) = A cos2#. 



(17) 



Below we will calculate the Josephson critical cur- 
rent Ic analytically and discuss the possibility of 
the atomic scale 0-ir transition. 

4. Josephson critical current 

The Josephson critical current Ic for a c-axis d- 
wave junction through the ferromagnetic insulator 
[Fig. 1(b)] can be expressed as 

Ic ^%J2J2 *K fe ' k 'M k > k ') F ( k > w„)F(fe', w„) 
fe,fe' »« 

(18) 

by assuming Ai = A 2 = A cos2# and thus F\ — 
F2 = F. We also assume that the tunneling matrix 
element t a (k, k') oc T a is given in terms of the co- 
herent tunneling in which the momentum fen paral- 
lel to the layer is conserved [22,23] and has the same 
Lp dependence on the transmission coefficient T a as 
Eq. (1) [7],i.e., 



tl(k,k%(k,k') = \t \ 2 (-l) LF 6 kllk , l 



(19) 



we obtain an analytical expression of Ic for T = 
K as 



Ic = 



where 



G 



N 



n 
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(20) 



(21) 



is the normal conductance with N being the density 
of states at Ep. The sign of Ic becomes negative 
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for the odd number of Lp and positive for the even 
number of Lp as was numerically found in [7-12]. 
Therefore, the spin dependent 7r-phase shift of the 
tunneling matrix element t a in the FI barrier gives 
rise to the atomic scale 0-n transition. 

5. Summary 

To summarize, we have studied the Joscphson 
effect in S/FI/S junction by use of the tunneling 
Hamiltonian method. We have analytically calcu- 
lated the Josephson current and showed the possi- 
bility of the formation of the atomic scale 0-ir tran- 
sition in such systems. This observation is consis- 
tent with previous numerically results. We hop that 
such FI based 7r-junctions become an element in the 
architecture of quantum information devices. 
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